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ABSTRACT 

6 ■ 

■ We have derived the column densities of heavy elements in three gamma-ray 

burst (GRB) optical transients, associated with the circumburst or interstellar 
medium (ISM) of the host galaxy. In comparison with the same elements observed 
^ ■ in damped Lyman-a (DLA) systems along QSO sight lines, we find evidence for 

much higher column densities of Zn II. The gap between the QSO-DLA and 
GRB-DLA distributions is smoothly bridged by observations of the interstellar 
absorption in the Milky Way and the Magellanic Clouds. Very small [Fe/Zn], 
[Si/Zn], and [Cr/Zn] values in GRB-DLAs indicate large dust depletions. Once 
the dust-to-metals ratios are determined, we find an optical extinction Ay ~ 1 
mags, to be compared with typical Ay < 0.1 in most QSO-DLAs. Our inference 
of high dust content appears to be in contradiction with the typical low reddening 
previously found in GRBs. One possible way to reconcile is a dust grain-size dis- 
tribution biased towards big grains, which would give a grey extinction. Possibly 
the small dust grains have been destroyed by the GRBs themselves. Our find- 
ings support the idea that primarily optically selected QSOs probe mainly low 



On leave of absence from Osservatorio Astronomico di Roma, Italy 



- 2 - 



gas/dust regions of high redshift galaxies, while the more powerful GRBs can be 
detected through denser regions of their ISM (molecular clouds and star forming 
regions). Therefore GRB-DLAs and QSO-DLAs together provide a more com- 
plete picture of the global properties of the interstellar medium in high redshift 
galaxies. 

Subject headings: cosmology: observations - gamma rays: bursts - galaxies: 
abundances - ISM: dust, extinction 

1. INTRODUCTION 

Gamma-ray bursts (GRBs) are among the most dramatic events in the Universe. Their 
luminosities are so high (10 51 — 10 54 ergs are emitted in a few seconds) that they can be seen 
at the highest redshifts. However, since the decline in their energy emission is very fast, the 
investigation and understanding of their physical nature is is not particularly easy (Piran 
2001). For instance, the redshift of the afterglow or host galaxy has been measured in only 
24 cases 2 . 

Only eight of these objects have been targeted quickly enough by ground-based tele- 
scopes to obtain low resolution spectra of the optical transient, and all show a common 
peculiarity: very strong UV absorption lines of low ionization species (Fiore 2001; Fynbo et 
al. 2002). This is similar to what is seen in damped Lyman-a systems (DLAs) associated 
with intervening galaxies along QSO sight lines. GRB-DLAs are certainly an unexplored 
and important source of information on the circumburst and/or interstellar medium (ISM) 
of the GRB host galaxy. Here we explore in more detail their properties and compare them 
with those of QSO-DLA galaxies. 

We use the curve of growth (COG) analysis to derive column densities of low ionization 
species in several GRB-DLAs. As for QSO-DLAs, we neglect ionization effects to derive 
the abundances of several heavy elements. The low ionization of the gas is warranted by 
the large equivalent widths of the low ionization lines, and the low equivalent widths of 
the moderate and high ionization lines. The difference in abundances between weakly and 
greatly dust depleted elements provides an estimate of the dust content. 

The COG analysis was applied to QSO-DLAs (Blades et al. 1982; Turnshek et al. 1989) 
before the advent of high resolution spectroscopy, which allows a direct fitting of individual 
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absorption lines. Since GRB emission fades very rapidly, high resolution (FWHM ~ 10 
km s" 1 ) spectra of GRB afterglows are very hard to obtain. 

In this work, we interpret strong differences between GRB-DLAs and QSO-DLAs as 
the result of the difference in the densities of the intervening gas. Since GRBs are, at least for 
some time, brighter than primarily optically selected QSOs, they can be observed even if their 
sight lines cross very dense and dusty regions, such as molecular clouds, where extinction is 
very severe. On the other hand, QSO sight lines are easier to detect when intersecting less 
dusty interstellar clouds. 

This bias in QSO-DLAs has been already suggested in the past (see, for instance, Fall 
& Pei 1993; Boisse et al. 1998). More recently, Ellison et al. (2001) have studied a small 
sample of radio-selected QSOs and tentatively found more DLAs than in optically selected 
QSOs. Even if this result is affected by small number statistics (significant at ~ 2a level), 
it supports the idea of dust obscuration in optically selected QSO-DLAs. 

Our analysis is alternative to and independent of the dust content estimate obtained so 
far for GRBs. The dust reddening in GRB afterglows has been derived in several objects (see 
for instance Reichart 1998; Jensen et al. 2001; Rhoads & Fruchter 2001; Fynbo et al. 2001; 
Galama & Wijers 2001) and in most cases it has been found to be small or negligible. The 
method applied assumes initial power law emission from the GRBs, and steep extinction 
curves such as those in the Milky Way (MW) or Large & Small Magellanic Clouds (MC). 
Possible discrepancies with the results of our analysis challenge the two initial assumptions 
of very steep extinction and initial power law emission of the central source. 

For this study, we examine the absorption lines of three of the eight GRBs with UV 
rest frame spectra, namely GRB990123 (z = 1.601, Kulkarni et al. 1999), GRB010222 (z = 
1.475, Masetti et al. 2001; Jha et al. 2001; Salamanca et al. 2001; Mirabal et al. 2002), and 
GRB000926 (z = 2.0379, Castro et al. 2001). We exclude GRB970508 (z = 0.835, Metzger 
et al. 1997), GRB990510 (z = 1.619, Vreeswijk et al. 2001) GRB000301C (z = 2.04, Jensen 
et al. 2001), GRB990712 (z = 0.4331, Vreeswijk et al. 2001), and GRB011211 (Holland et 
al. 2002) for poor quality of the spectra and therefore lack of interesting absorption features. 

2. GRB-DLAs COLUMN DENSITY DETERMINATION 

Similarly to QSO-DLAs, GRB-DLAs show strong Fe II, Mg II, and Si II absorptions, 
together with relatively weak Zn II and Cr II. Table 1 lists oscillator strengths f\ (from 
the recent compilation by Prochaska et al. 2001) and rest frame equivalent widths W r of 
absorption lines in the three GRBs used in the present analysis. The large W r of all these 
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species relative to higher ionization species like C IV, Si III and Si IV, indicate that the 
medium is weakly ionized. From the equivalent widths (EWs), we derive column densities 
and directly compare the results with typical column densities in QSO-DLAs. Comparing 
column densities instead of metallicities [X/H] has the advantage that the H I column density 
(known only in GRB000926 with large uncertainty; Fynbo et al. 2002) is not required. 
Moreover, it is not necessary to make any assumptions about the dust depletion, cosmic 
abundance or ionization level of the considered element, but we can still search for similarities 
and differences between GRB-DLAs and QSO-DLAs. 

A direct comparison is complicated by the fact that normally QSO-DLAs are observed at 
high spectral resolution (FWHM ~ 10 km s _1 ) in a large wavelength range (thanks to echelle 
spectroscopy), allowing a good determination of column densities of many ions, even in cases 
where saturated or weak lines in complex features are detected (see for instance Lu et al. 1996; 
Pettini et al. 1997; Prochaska & Wolfe 1999). On the other hand, the typical resolution of 
GRB optical spectra is FWHM > 100 km s -1 . Moreover, the majority of absorption lines 
are saturated, and it is not possible to use the linear part of the COG to determine column 
densities. However, the lines are numerous and so, even if the uncertainty may be large, it 
is possible to make an abundance analysis using the general COG analysis (Spitzer 1978). 
Since the typical resolution of GRB spectra does not allow the separation of the different 
components of a complex feature, due to the many clouds along the line of sight, in our 
analysis we adopt a large "effective" Doppler parameter, being the result of the superposition 
of many narrow absorption lines. For reference, we note that typical velocities in QSO- 
DLAs are ~ 100 km s -1 or more. The COG analysis applied with the single component 
approximation was used for QSO-DLA studies before high resolution spectroscopy of QSOs 
became possible (see for instance Blades et al. 1982; Turnshek et al. 1989). In fact, Jenkins 
(1986) showed that the COG technique applied to complex features gives nearly the correct 
answer (the simulated-to-true column density ratio rarely goes below 0.8) even if different 
lines have very different saturation levels or Doppler parameters. This general result is also 
confirmed by Savage & Sembach (1991). 

The best-fitting column density is determined by using simultaneously different EWs 
of various ions and choosing the COG with an effective Doppler parameter that minimizes 
the x 2 - Once the best effective Doppler parameter is determined, the error on the column 
density for each ion is determined by varying log A" for which x 2 /dof < 1. In the following 
§ 2.1-2.3, we discuss column density measurements in each individual GRB. Figures 1-3 
show the best -fitting W r /X and column densities, and the corresponding COG. Estimated 
column densities are given in Table 2. 
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2.1. GRB990123 

The optical transient of GRB990123 at z = 1.6004 was observed with the Keck/LRIS 
(Kulkarni et al. 1999) for a final resolution and signal-to-noise ratio of FWHM = 11.6 
A (400 - 700 km s" 1 ) and S/N « 30, respectively. The very weak Fe II A2260 EW gives 
inconsistent results with other Fe II lines, and thus this detection has not been considered in 
our analysis. Fe II A2374 is not reliable because it is probably partly blended with the Fe II 
A2382 line. 

For the four remaining Fe II lines, we get a best -fit of logiV(Fe II) = 14.78^q;}q and 
b = 50 km s" 1 (Figure 1). From the EW of Zn II A2026 and considering b = 50 km s _1 , we 
derive a Zn II column density of log iV(Zn II) ~ 14, and calculate a Zn II A2062 contribution 
to the Zn II A2062 + Cr II A2062 blend of W r = 0.58 A, to be compared with the detected 
W r = 0.48 ±0.04 A. This means that the contribution of the Cr II A2062 line to the blend is 
negligible or non existing. Therefore we recalculate the Zn II column density from the Zn II 
A2026 and Zn II A2062 doublet, from which we derive log iV(Zn II) = 13.95+°;^ 



2.2. GRB000926 

GRB000926 at z = 2.0379 was observed with the Keck/ESI echelle spectrograph (Castro 
et al. 2001), with a good signal-to-noise ratio (S/N = 10 — 20) and relatively high resolution 
(FWHM ~ 80 km s _1 ). Castro et al. (2001) also show a low resolution spectrum, with EW 
measurements that are systematically larger than those obtained from the high resolution 
spectrum, but the authors claim that the high resolution measurements are more accurate, 
so we only use these. The spectrum shows at least two absorbing clouds separated by 168 
km s -1 . In our analysis, we consider the total equivalent widths of the two components 
together, since the resolution does not allow them to be clearly separated. 

The observed Si II A1808 and Si II A1526 lines, and the five detected Fe II lines give 
logiV(Si II) = 16.47+^ and logiV(Fe II) = 15.6±g;?g, respectively (b = 115 km s" 1 , Figure 
2). If we assume b = 115 km s" 1 , we get logiV(Cr II) = 14.34+°:os from the Cr II A2056 
EW. This corresponds to W r = 0.5 - 0.6 A and W r = 0.35 - 0.43 A for Cr II A2062 and 
Cr II A2066, respectively. Neither of these measurements is consistent with those observed 
(W r < 0.28 A and W r < 0.2 A, respectively); this might indicate a contamination of the 
Cr II A2056 line by an unidentified absorption associated with an intervening metal system. 
Moreover, the Cr II A2062 is probably corrupted by a noise spike in the spectrum (see Figure 
3 of Castro et al. 2001) probably due to a remnant of the strong sky line at A = 6258.05 
A. The same contamination might also effect the Zn II A2062 line (the EW is inconsistently 
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lower than expected from Zn II A2026), therefore the Zn column density is determined from 
Zn II A2026 only: logiV(Zn II) = 13.82 ± 0.05 for b = 115 km s _1 . The lack of Fe II A2586 
detection is puzzling. If, for instance, we assume W r (Fe IIA2586) < 1 A (way above the 
detection limit in the spectrum), the line is almost unsaturated, and log iV(Fe II) < 14.6 for 
b > 70 km s _1 , inconsistent with what is found from the other Fe II lines. 

Fynbo et al. (2001) report a tentative H I column density measurement of 7V(H I) 
2 x 10 21 cm -2 . This leads to a relatively high metallicity with [Zn/H] ~ —0.13. 

2.3. GRB010222 

The optical transient of GRB010222 at z = 1.475 was observed with the FLWO 1.5m 
telescope (Jha et al. 2001; FWHM = 6 A or 300 - 450 km s -1 , S/N ~ 10), with the 3.58m 
TNG (Masetti et al. 2001; FWHM = 4.8 A or 200 - 400 km s _1 , S/N = 10 - 20), with 
the 4.2m WHT (Salamanca et al. 2001; FWHM = 3.3 - 5.8 A or ~ 300 km s _1 , S/N ~ 10), 
and with Keck/LRIS and ESI (Mirabal et al. 2002; FWHM = 11 - 13 A or ~ 650 km s~\ 
and FWHM = 0.4 — 0.8 A or ~ 30 km s _1 ). Even though the different spectra span a 
time interval of ~ 27 hours, starting 5 hours after the burst, the absorption line equivalent 
widths do not show significant time variability (Mirabal et al. 2002). Therefore, to obtain a 
better estimate of column densities, we combined the EWs of the same lines from different 
observations weighted according to the errors. We constrain the effective Doppler parameter 
to b = 70 km s^ 1 , using the Si II A1526 and Si II A1808 lines, and seven Fe II lines (Figure 
3). We get logiV(Si II) = 16.09 ± 0.05 and logiV(Fe II) = 15.32°;^. The Fe II lines are 
scattered in the COG diagram, probably because of non-uniform EW measurements by the 
four groups. 

The Zn II A2026 line at A ss 5018 A is blended with the Mg I A2026 line. The contam- 
ination of Mg I to this doublet is estimated using the Mg I A2852 line at A ~ 7065 A and 
W r (Mg I A2852) = 1.22 ± 0.04 A. Assuming an effective Doppler parameter b > 35 km s _1 , 
we get logiV(Mg I) < 14 and W r {Mg I A2026) < 0.04 A. We then assume that the feature 
at A « 5018 A is dominated by the Zn II A2026 absorption, with W r = 0.79t f 7 A and the 
lower error determined by a possible contamination from Mg I A2026. 

The Cr II A2056 EW is very low, therefore it must be weakly or not at all saturated. 
From Cr II A2056 and Cr II A2066 lines, we get logiV(Cr II) = 14.04+°;^. From this, we 
calculate a contribution of Cr II to the Zn II A2062 + Cr II A2062 doublet of W r (Cr II 
A2062) = 0.32 ± 0.04 A. For this doublet we only use the measurements given by Jha et 
al., Salamanca et al. and Mirabal et al. (weighted mean W r = 0.73 ± 0.02 A), because the 
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spectrum by Masetti et al. shows strong blending with other absorption lines. Once the 
Cr II contribution is taken into account, we derived W r (Zn II A2062) = 0.41 ±0.05 A. Using 
also the Zn II A2026 detection, we obtain logiV(Zn II) = 13.78 ± 0.07 for b = 70 km s _1 . 
The Mn II column density is obtained from the Mn II AA2576, 2594, 2606 triplet. For b = 70 
km s" 1 , we get log N{Mn II) = 13.61g;g|. 



3. HEAVY ELEMENT ABUNDANCES IN GRB-DLAs 

Figure 4 displays the QSO-DLA column density histograms 3 of Fe II, Cr II, Zn II, and 
Si II. These column densities are obtained from a large compilation of QSO-DLA measure- 
ments collected from the literature (see Savaglio 2000 for a description of the compilation). 
The sample contains 98 QSO absorption line systems associated with mostly neutral gas 
clouds (H I column density larger than ~ 10 19 cm~ 2 ) of the ISM in 0.0 < z < 4.6 galaxies, 
for which the column density of one or more of the following elements is measured: Zn II, 
Si II, Cr II, Fe II. The symbols in Figure 4 indicates the column densities of the same ions 
measured for the three GRB-DLAs. This reveals several differences between QSO-DLAs and 
GRB-DLAs. The Zn II column density in GRB-DLAs shows the most striking difference, 
being much larger than in QSO-DLAs. For Si and Cr the deviations are also very large. 
Fe II behaves differently, being quite consistent with the upper part of the Fe II QSO-DLA 
distribution. 

If the gas is nearly neutral, as suggested by the large EWs of the low ionization lines, 
the ionization correction can be neglected and it is possible to calculate the element relative 
abundance from ions with ionization potential above 13.6 eV. This assumption is also widely 
used to derive abundances in QSO-DLAs. Table 3 shows the comparison between heavy 
element relative abundances in QSO-DLAs (mean values, column four) and in the three 
GRB-DLAs (columns five, six and seven) . Columns two and three also give the number of 
QSO-DLAs considered and the mean redshift, respectively. For both GRB-DLAs and QSO- 
DLAs, the standard solar abundances from Grevesse, Noels, & Sauval (1996) are adopted. 
The very low [Fe/Zn] values in three GRB-DLAs (-2.03, -1.08, and -1.32) are much smaller 
than in QSO-DLAs: ([Fe/Zn]) = —0.46 ±0.24 in 32 objects. This is an indication of strong 
dust content in the GRB circumburst medium and/or in the ISM of the host galaxy, because 
while Fe is heavily depleted on dust in the ISM of galaxies, Zn is marginally depleted only in 



3 We have not used Al II measurements detected in GRB000926 and GRB010222, because not included 
in the QSO-DLA sample, and because [Fe/Al] = -0.3 ±0.3 and -0.4 ±0.2 for GRB000926 and GRB010222, 
respectively, which is not particularly interesting for our discussion. 
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the densest clouds (Savage & Sembach 1996a). The plot of [X/Zn] abundances as a function 
of the Zn II column density is also very interesting (Figure 5). In QSO-DLAs these show a 
trend, the [X/Zn] being smaller for larger log iV(Zn II). The [X/Zn] values for GRB-DLAs are 
located at the extension of the [X/Zn] vs. log iV(Zn II) distribution in QSO-DLAs. The gap 
between QSO-DLAs and GRB-DLAs is filled by a completely different class of absorbers: 
those in the ISM of the MW and MC. Low values of [Fe/Zn] are also an indication of high 
dust content. Petitjean, Srianand, & Ledoux (2002) discuss the gas density in QSO-DLAs 
with H 2 detection, and suggest that the dust depletion might be larger in denser clouds. If 
this can be extended to GRB-DLAs, it might imply that their gas density is higher than in 
QSO-DLAs. 

The small [X/Zn] values found in GRB-DLAs might be a general property of high 
metal column density systems rather than a peculiarity of GRBs themselves. GRBs are 
initially brighter than QSOs, hence detectable even in cases of heavy dust obscuration. In 
other words, QSO-DLAs and GRB-DLAs might probe the same population of galaxies, but 
different regions: the sight lines of the former cross moderately dense regions, while the 
sight lines of the latter might cross denser regions such as molecular clouds in a star forming 
environment. 

As already discussed in § 2, the COG analysis may underestimate column densities if 
the lines are strongly saturated, as in the case of Fe II. This would artificially create a large 
dust depletion if [Fe/Zn] were used as an indicator. However, the underestimate of saturated 
lines is in most cases not larger than 0.1 dex (Jenkins 1986). Moreover, an independent check 
of consistency is provided by the comparison of Fe II with Si II and Cr II. Si II and Fe II 
have similar problems with saturation, but Si is much less dust depleted (Savage & Sembach 
1996a). The iron-to-silicon relative abundance is low in the two GRB-DLAs discussed here, 
indicating high depletion of iron in dust: [Fe/Si] = — 0.83±g;f| and -0.73^ in GRB000926 
and GRBO 10222, respectively. Also [Cr/Zn] is a dust indicator. Cr II absorption lines are 
generally weakly saturated, but the column densities in the same two GRB-DLAs are low 
compared to Zn II: [Cr/Zn] = -0.51 ± 0.07 and -0.77^;^ in the same two GRBs. 

There are indications that the H I column density in GRB-DLAs is larger than in typical 
QSO-DLAs. In 89 QSO-DLAs with log N(R I) > 20.0, only 15 (17%) have \ogN(R I) > 
21.0. The H I column density has been estimated in GRB000301C by Jensen et al. (2001), 
who report iV(H I) = 1.5^^ x 10 21 cm -2 . GRB000926 has also been observed with the Nordic 
Optical Telescope (Fynbo et al. 2002) and a very strong Lja absorption line identified, with 
column density ~ 2 x 10 21 cm~ 2 . This is in a very noisy part of the spectrum, but if we assume 
that this H I column density is close to the correct value and neglect ionization correction, we 
derive [Si/H] = -0.38, [Zn/H] = -0.13, [Cr/H] = -0.64, [Fe/H] = -1.21, suggesting a large 
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metallicity and dust depletion for a large H I GRB-DLA. In comparison, in 8 QSO-DLAs 
with logiV(H I) > 21.0, and measured Zn II and Cr II, we find ([Zn/H]) = -1.46 ± 0.34 and 
([Cr/H]) = -1.74 ±0.22. 



4. DUST DEPLETION AND EXTINCTION 



Since four elements are measured in GRB000926, it is possible to derive its dust depletion 
pattern. We adopt the method described in Savaglio (2000). Basically, we consider the four 
depletion patterns observed in the Milky Way, e.g. the depletion patterns in the warm halo 
(WH), warm disk+halo (WDH), warm disk (WD), and cool disk (WC) clouds (Savage & 
Sembach 1996a). For each depletion pattern, we find the best-fit values for the dust-to- 
metals ratio (d/dj, where j depends on the depletion pattern assumed) and the metallicity 
compared to solar Z/Z & . The result for GRB000926 at z — 2.038 is shown in Figure 6. The 
best fit is given by the warm halo cloud pattern, with Z/Z & = 0.78. This is higher than the 
metallicities in QSO-DLAs: (Z/Z Q ) = 0.22 ± 0.23 in 21 DLAs at 1.7 < z < 2.3 (Savaglio 
2000). 

We can now estimate the obscuration of the GRB due to dust distributed along its line 
of sight. This is proportional to the column density of metals. The extinction in the V-band, 
due to dust, can be approximated by: 



nr d z N(m) 

A V = 0.5A- — v J _ 2 1 

dwH Z e 10 21 cm 1 

where 0.54 is the typical value of Ay in the solar neighborhood for a column density of gas 
with log JV(H I) = 21.0. For GRB000926, we find A v « 0.9 mag. This is much higher than 
in QSO-DLAs, where V extinction is typically lower than 0.1 mag. 

The depletion pattern and extinction can also be calculated for GRB990123 and GRB010222, 
because this depends on the product of metallicity and H I column density, i.e., the column 
density of metals. We consider the Zn II column density and neglect the dust depletion 
correction for this element. Best-fits for these two GRB-DLAs give a warm halo cloud de- 
pletion pattern (Figure 7), and dust extinctions of A v 1.1 and 0.7 mags for GRB990123 
and GRB010222, respectively. 
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5. DISCUSSION 

We have used the curve of growth to study the heavy element absorption lines asso- 
ciated with three GRBs, namely GRB990123 (Kulkarni et al. 1999), GRB000926 (Castro 
et al. 2001), and GRB010222 (Masetti et al. 2001; Jha et al. 2001; Salamanca et al. 2001; 
Mirabal et al. 2002). These are the only GRBs with available UV spectra obtained so far 
with quality good enough to detect high S/N absorption lines. All GRB spectra show very 
high Zn II absorption associated with the circumburst medium or interstellar medium of 
the host galaxy, from which we derive a mean column density of (logiV(Zn II)) = 13.85 at 
(z) = 1.70. This is ~ 3 times larger than the largest Zn II column density detected so far in 
QSO-DLAs. High Zn II columns indicate high metallicity and/or low ionization of the gas, 
also supported by the presence of strong Fe II, Si II, and Cr II, as well as Mg II together 
with Mg I. Moreover, the GRB-DLAs analyzed here satisfied the criteria found by Rao & 
Turnshek (2000) that are used for QSO-DLA selection, namely W r (Mg II A2796) > 0.6 
A and W r (Fe II A2600) > 0.5 A. 

Since the metal column densities are so high, the dust extinction can be very important. 
A good indicator of dust content is given by the [Fe/Zn] values; in the GRB-DLAs these are 
much smaller (higher dust depletion) than those seen in QSO-DLAs: ([Fe/Zn]) = —1.48 in 
the former, ([Fe/Zn]) = —0.46 in the latter. The combination of large column densities of 
metals and large dust depletion can lead to large extinction in GRB-DLAs. If we compare 
metal column densities with the ones in the solar neighborhood, where extinction is very well 
studied, we find optical extinction of Ay m 0.7, 0.9 and 1.1 magnitudes. This is much larger 
than that found for most QSO-DLAs: Ay < 0.1. We notice, however, that another possible 
explanation for the very low [Fe/Zn] values would be a nucleosynthetic production of Fe and 
Zn in GRBs and/or their environment being different from the standard solar pattern. 

These strong differences between QSO-DLAs and GRB-DLAs do not necessarily indi- 
cate a difference in the nature of the host galaxies. Indeed, the distribution of the Zn relative 
abundances as a function of Zn II column density in GRB-DLAs and QSO-DLAs is consis- 
tent for these two classes of objects. [Fe/Zn], [Cr/Zn], and [Si/Zn] decrease as a function 
of logiV(Zn II) for QSO-DLAs and GRB-DLAs, which are located in the upper left part 
and in the lower right part of the diagram, respectively. The QSO-DLA and GRB-DLA 
distributions are bridged by the heavy element abundances in the diffuse and dense ISM 
of the Milky Way and the Magellanic Clouds. This suggests that the galaxies associated 
with the absorbing gas around GRBs are not different from those of QSO-DLAs. The im- 
portant difference is in the column densities of the intervening ISM: GRBs are (at least for 
some time) much brighter than optically selected QSOs, so their sight lines can probe much 
denser and dustier regions of galaxies. Indications of a bias effect in QSO-DLAs due to dust 
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obscuration have been pointed out by several authors (Fall & Pei 1993; Boisse et al. 1998; 
Savaglio 2000; Petitjean et al. 2002). 

Reichart (1998), Jensen et al. (2001), Rhoads & Fruchter (2001), Fynbo et al. (2001), and 
Galama & Wijers (2001) have estimated the dust extinction in GRBs from curvature in their 
spectra, assuming that the intrinsic GRB spectra from the UV to IR bands are power laws. 
This approach gives a small V extinction only if MW or MC laws are assumed: Ay m — 0.2 
for GRB970508, GRB000301C, GRB000926, and GRB990123. The result for GRB990123 
is supported by the lack of spectral slope variation going from the optical to x-ray bands 
(Castro-Tirado et al. 1999). Simon et al. (2001) have found that the color variation during 
GRB decline in 17 objects (three of which are the ones studied in this paper) is relatively 
small, suggesting similar and low reddening in these GRBs. However, they also notice that 
another sample of five GRBs (not included in the first one, for lack of complete photometric 
information) have much redder colors. Low extinction can be explained by sublimation of 
dust grains due to UV emission (Waxman & Draine 2000). Fruchter, Krolik, & Rhoads 
(2001) suggest grain destruction by GRB heating and grain charging. These effects can be 
present up to a distance of ~ 100 pc from the GRB source. 

This is in principle in contradiction with our findings. However, the dust disruption 
mechanisms suggested by Fruchter et al. (2001), or other mechanisms, such as sputtering 
(Draine & Salpeter 1979), predict that small grains are destroyed first. Dust destruction of 
small grains would still give a large content of dust, if the initial grain-size distribution is 
dominated by large grains, as suggested in models of MW an MC extinction curves (Wein- 
gartner & Draine 2001). Therefore, what remains in dust can still be an important fraction 
of the total metal content. Grain coagulation favored by high densities, leading to the for- 
mation of large grains, is also possible (Mathis 1990). In these cases, MW or MC extinction 
curves might be inappropriate to describe dust obscuration. A dust grain distribution skewed 
toward large grains would give a grey extinction, therefore may explain the low reddening 
found in GRB afterglows (Stratta et al. 2002). This is a plausible speculation, but quanti- 
tative calculations, needed to determine if it works in practice, are beyond the scope of this 
paper. A first attempt is instead presented by Perna & Lazzati (2002), who have developed 
a detailed numerical code that describes, among other effects, the dust destruction under 
an intense X-ray and UV radiation field the very first seconds after the burst. Alternative 
explanations are unknown systematic effects. For instance, the low reddening suggested re- 
lies on the assumption of an intrinsic power law emission of the optical afterglow, which in 
principle cannot be tested. It might also be that the extinction in itself is a power law, and 
when combined with the intrinsic GRB power law emission, would still give a power law 
spectrum. Therefore, this would require low or no rectification. 
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The main result of the present work is the strong difference between QSO-DLA and 
GRB-DLA heavy element column densities. GRB-DLAs occur most likely along the sight 
line of much brighter sources, and are embedded in star forming and therefore metal rich 
regions, where the dust extinction is probably not negligible. This is another evidence 
supporting the idea that DLAs in front of the fainter QSO radiation, can be detected and 
studied in detail only when the sight line crosses less dense, metal polluted and/or extincted 
regions of high redshift galaxies. 

We thank Tim Heckman, Fiona Harrison, Eliana Palazzi and Ken Sembach for fruitful 
information on GRB spectra and ISM absorptions. We are particularly grateful to Daniela 
Calzetti, Julian Krolik, Nicola Masetti, Nino Panagia, and James Rhoads for stimulating 
discussions that helped to improve this paper. 
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Fig. 1. — Curve of growth for GRB990123 absorption lines (b is the effective Doppler param- 
eter in km s _1 ). Points give the best-fit column densities, reported in Table 2. The straight 
line is the linear case (b = +oo). 
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Fig. 4. — Histogram column densities of Fe II, Cr II, Zn II, and Si II in QSO-DLAs. The 
GRB-DLA column densities are reported as filled triangles (GRB990123), empty squares 
(GRB010222) and filled circles (GRB000926). 
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Fig. 5. — Iron-to-zinc, chromium-to-zinc, and silicon-to-zinc relative abundances, vs. Zn II 
column density, in QSO-DLAs (empty squares), in the ISM of the Milky Way and the 
Magellanic Clouds (crosses), and in GRB-DLAs (filled symbols). Data for the local ISM 
are taken from Cardelli, Sembach, & Savage (1995), Roth & Blades (1995, 1997), Savage 
& Sembach (1996b), Sembach & Savage (1996), Spitzer & Fitzpatrick (1995), and Welty et 
al. (1997, 1999). Mean errors for [X/Zn] in QSO-DLAs is ~ 0.08 dex. Errors for local ISM 
column densities, when available, are typically < 0.1 dex. 
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Fig. 6. — Depletion patterns in the absorbing gas of GRB000926. The models (lines) are 
taken from average gas-phase abundance measurements in warm halo, warm disk+halo, 
warm disk, and cool disk clouds of the Milky Way (Savage & Sembach 1996a). The metal- 
licity, dust-to-metals ratio, and best-fit x 2 are also given. 
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Fig. 7.— Best-fit depletion patterns in the absorbing gas of GRB990123 and GRBO 10222. 
The models (lines) are taken from average gas-phase abundance measurements in warm disk 
clouds of the Milky Way (Savage & Sembach 1996a). Dust-to-metals ratio and best-fit x 2 
are also given. 
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Table 1: Rest-frame equivalent widths 
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a Cr II A2062 contribution to the Zn II A2062 + Cr IIA2062 blend is negligible. 

b Noisy region of the spectrum, given only upper limit based on the total EW + 3er, reported by Castro ct 
al. 2001. 

c Weighted mean from Jha et al. 2001, Salamanca et al. 2001, and Mirabal et al. 2002. 
d Lower error found assuming possible contamination from Mg I A2026 line. 

Contribution of the Cr II A2062 line to the Zn II A2062 + Cr IIA2062 blend calculated using the Cr II 
column density as derived from Cr II AA2056, 2066. 
•^Weighted mean from Masetti et al. 2001, Salamanca et al. 2001, Jha et al. 2001, and Mirabal et al. 2002. 
^Weighted mean from Salamanca et al. 2001 and Mirabal et al. 2002. 
^Mirabal et al. 2002. 

l Weighted mean from Masetti et al. 2001, Jha et al. 2001, and Mirabal et al. 2002. 
•'Weighted mean from Jha et al. 2001 and Mirabal et al. 2002. 
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Table 2: Column densities 
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Table 3: Element relative abundances in QSO-DLAs and GRB-DLAs 
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